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Abstract 

We present results of our numerical calculation of the mass spectrum 
for isospin one-half and spin one-half non-strange baryons, i.e. the ground 
and excited states of the nucleon, in quenched lattice QCD. We use a new 
lattice discretization scheme for fermions, domain wall fermions, which pos- 
sess almost exact chiral symmetry at non-zero lattice spacing. We make a 
systematic investigation of the negative-parity A^* spectrum by using two dis- 
tinct interpolating operators at /3 = = 6.0 on a 16^ X 32 X 16 lattice. 
The mass estimates extracted from the two operators are consistent with each 
other. The observed large mass splitting between this state, A^*(1535), and 
the positive-parity ground state, the nucleon A^(939), is well reproduced by 
our calculations. We have also calculated the mass of the first positive-parity 
excited state and found that it is heavier than the negative-parity excited 
state for the quark masses studied. 
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I. INTRODUCTION 



An important challenge in lattice calculations is to reproduce the hadron mass spec- 
trum from first principles in quantum chromodynamics (QCD). The latest lattice QCD 
calculations of the light-hadron mass spectrum in the quenched approximation agree with 
experimental values within about 5% |T]-§|. However, this success is mainly restricted to 
ground states. Indeed, results are scarcely available for the excited-state mass spectrum. 

Another essential shortcoming of these calculations which use Wilson or Kogut-Susskind 
fermions is the absence of chiral symmetry at finite lattice spacing, in accord with the 
Nielsen-Ninomiya no-go theorem ^ . At non-zero lattice spacing Wilson fermions explicitly 
break the full chiral symmetry of the continuum down to the vector sub-group, so only 
flavor symmetry is preserved. On the other hand, Kogut-Susskind fermions have only a 
single exact U{1) axial symmetry, and flavor symmetry is completely broken Of course, 
it is expected that in the continuum limit, which is difficult to achieve in practice, both 
actions recover the full chiral symmetry. 

Several years ago Kaplan constructed a new type of lattice fermion known as domain 
wall fermions, which were further developed by Shamir and also by Narayanan and 
Neuberger 0. Especially, the former reformulated it for lattice QCD simulations. The key 
feature of domain wall fermions is that they utilize an extra fifth dimension to circumvent the 
Nielsen-Ninomiya no-go theorem and maintain chiral symmetry at non-zero lattice spacing. 
In practical simulations the extra dimension is finite, so the chiral symmetry is not exact. 
The symmetry breaking is very soft, however, since it is highly suppressed with the number 
of sites in the extra dimension, Lg. In other words, Lg gives us a way to control the violation 
of chiral symmetry. Domain wall fermions also possess exact flavor symmetry for any value 
of Lg. 

Quenched lattice QCD calculations with domain wall fermions have shown that good 
chiral properties are obtained for moderate sizes of the fifth dimension, ~ 10 — 16, if the 
lattice spacing is small enough (a ^ 0.1 fm) [[TO| , p!T| . Recent studies by the RIKEN-BNL- 
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Columbia-KEK collaboration [0] and the CP-PACS collaboration |T3| have quantified in 
detail the explicit chiral symmetry breaking effects due to finite L^. For low energy QCD, the 
results can be simply summarized: there is a unique additive quark mass, mres, appearing in 
the low energy QCD lagrangian which arises from the finite size of the extra dimension [jl2[ . 
This additive quark mass has been measured quite accurately, and for Lg = 16 and quenched 
f3 = 6.0, mres ~ 3% of the strange quark mass [l^JT^. In this paper, we apply domain wall 
fermions to baryon excited states, especially the spin one-half and isospin one-half negative- 
parity nucleon, A^*(1535), as a further test of domain wall fermions in the baryon sector. For 
masses which are 0(1) GeV in the chiral limit, we do not expect mres to have a significant 
effect. 

We are interested in a long-standing puzzle in the excited state spectrum of the nucleon. 
The first question addressed in this paper is whether the mass difference between the nucleon 
A^(939) and the negative-parity nucleon A^*(1535) is well reproduced in lattice QCD. The 
spin one-half A^* state can be considered the parity partner of the nucleon. Of particular 
interest is the large mass splitting between and A^*. From the viewpoint of parity partners, 
these two states would be degenerate if the relevant 2-flavor chiral symmetry were exact and 



preserved by the vacuum WM. Of course, there is no proof that this large mass splitting comes 
directly from the spontaneous breaking of chiral symmetry. However, at least spontaneous 
chiral symmetry breaking is responsible for the absence of such parity doubling since the 
explicit breaking is quite small in the case of two flavors. In this sense, regardless of a model 
or a theory, chiral symmetry and its spontaneous breaking are important for reproducing 
precisely the mass splitting between parity-partner hadrons. In this paper we show that 
domain wall fermions accurately reproduce the large observed mass splitting (some of our 
results have been reported earlier |]15|,|16[). 

Conventional lattice fermion schemes have had difficulty in this challenge. An early 



calculation using Wilson fermions as well as a more recent one |T8[, both just managed 
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to extract a mass splitting between the parity partners with large uncertainties.^ 



The first sophisticated calculation []T9[ using an improved Wilson fermion action and 
relatively heavy quarks resolved a definite mass splitting between the parity partners which, 
however, was about a factor of two smaller than the observed splitting. Recent results 
using improved Wilson fermions have confirmed the large mass splitting over a wide range 
of quark masses that we have found with domain wall fermions pO| , ^ . We note that 
the leading lattice spacing errors that are removed from these improved calculations break 
chiral symmetry. Although Kogut-Susskind fermions have a remnant U{1) axial symmetry, 
they cannot be used practically for the spin one-half A^* mass calculation due to fiavor 
mixing. The reason is that Kogut-Susskind fermions have only discrete fiavor symmetries 
belonging to a subgroup of SU (4) which contains three irreducible representations, 8, 8' 
and 16 for baryon operators. Two appropriate representations 8 and 16, to which A^*(1535) 
belongs, also contain the negative-parity A states A(1405) and A(1520) and the spin three- 
half negative-parity nucleon state A^*(1520) 0. Thus, the study of the spin one-half A^* 
spectrum with Kogut-Susskind fermions always faces inevitable contamination from lower 
mass states. 

It is also interesting to note that a non-relativistic quark model with the so-called color 
magnetic interaction and the MIT bag model [^], both of which explicitly break chiral 



symmetry, have some difficulty reproducing the large mass splitting between A^(939) and 
A^*(1535) without adoption of less realistic model parameters. The non-relativistic quark 
model is based on a harmonic oscillator description of the orbital motion of constituent 



quarks. As remarked in the appendix of the original paper by Isgur and Karl , the plau- 
sible value of its oscillator quantum should be roughly 250 MeV to reproduce the observed 
charge radius and magnetic moment of the nucleon. Since this model regards A^*(1535) as 



^In a recent preprint [40| which appeared after this work was complete it is clear that the chiral 
symmetry of their fermions was significantly improved by simulating closer to the continuum limit 
than previous studies which used Wilson fermions and by adding the clover term. 
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a state with one quantum excitation in orbital motion, it indicates that the corresponding 
A^* state hes at most a few hundred MeV above the ground state. Even worse, there is the 
serious problem of the wrong ordering between A^*(1535) and the positive-parity excited nu- 
cleon A^'(1440) because the corresponding A^' state should be assigned two oscillator quanta 
in this model P2[. This wrong ordering problem does not seem to be easily alleviated 



It is easy to see that the MIT bag model faces essentially the same problem, as the single 
quark states in the model alternate in parity with roughly even spacings |^ . 

The question arises how does this ordering of A^* and A^' appear in lattice QCD calcu- 
lations? Until |jl5[ this question could not be answered for lack of systematic results and 
investigations. A calculation of the mass of the positive-parity excited nucleon is, of course, 
much more difficult than the nucleon ground-state. Attempts have been made to evaluate 



the A^' mass from a two state fit to the nucleon correlation function |25|. However, large 
statistics are required compared to a single exponential fit. Also it is difficult to control 
the systematic errors. Leinweber used the QCD sum rule (QCDSR) continuum model and 



quenched lattice data to estimate the threshold of contributions from excited states [26 



However, the scaling of his result has been questioned by Allton and Capitani [B^ . In this 
paper we take an alternative approach, using the continuum- like behavior of the domain 



wall fermion operator [|T5l to obtain the mass of the positive-parity excited nucleon. Our 



results have been confirmed in pO[ where a similar approach is employed in the context of 
improved Wilson fermions. 

In conventional lattice QCD calculations an interpolating operator which is strongly as- 
sociated with the non-relativistic limit is used to extract the mass of the nucleon ground 
state. A second unconventional operator, which does not have a non-relativistic limit, is dis- 
carded since it is expected to couple weakly to the ground state. Indeed, using this operator, 
no one has succeeded in evaluating the mass of the nucleon ground state in lattice QCD 
calculations with Wilson fermions I^H] because of its small coupling to the ground state |^ . 
The expectation of an approximately zero overlap on the ground state provides the possi- 
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bility that the use of the second operator in lattice calculations directly yields the mass of 
the A^' state, at least in the relatively heavy (valence) quark mass region. This prospect, 
however, is built on the assumption that the lattice defined operators inherit the features 
of the continuum ones. In the case of the Wilson fermions, the Wilson term, which explic- 
itly breaks chiral symmetry, induces mixing between the conventional and unconventional 



operators [^,^. Thus, the desired feature of the unconventional operator in the contin- 
uum is diminished in lattice calculations with Wilson fermions. On the other hand, this 
type of mixing between three quark operators is absent at one loop in perturbation theory 



using domain wall fermions with large Lg thus we expect this mixing to be suppressed 



in domain wall fermion lattice calculations. Indeed, we find that the second operator has 
negligible overlap with the ground state and furthermore provides a reliable signal of the 
excited state as the asymptotic state in the heavy quark mass region. 

The organization of our paper is as follows. In Sec. II, we briefly review the basic formulae 
and notation regarding domain wall fermions. In Sec. Ill we investigate the properties of 
the two-point correlation function for the nucleon and its parity partner. Sec. IV gives the 
details of our Monte Carlo simulations and the results for the parity partner of the nucleon 
(A^*) and the first positive-parity excited nucleon (A^'). Then, we compare our results to the 
experimental values. Finally, we present our conclusions in Sec. V. 

II. DOMAIN WALL FERMIONS 

In this section we closely follow the development of domain wall fermions by 
Shamir The domain wall fermion action is essentially regarded as a five- dimensional 

extension of the Wilson fermion action: 

Sdwf = - E E ^(^' ■^)[ks'Dl,, + 6,,,,D^J^{x', s') , (1) 

where x, x' are four-dimensional Euclidean space-time coordinates and s, s' denote coor- 
dinates in the extra dimension labeled from to — 1 (to take advantage of existing 
high-speed computer code, our domain wall fermion Dirac operator is the Hermitian conju- 
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gate of the one in hence our notation is the same as in 10). Here, dI^, corresponds to 
the four-dimensional Wilson-Dirac operator with a mass term (domain wall height) M5. 

Dly = ^ E {(1 - 7M)f^M(^) W' + (1 + 7M)^i(^')5x-A,x'} + (M5 - 4)5.,., , (2) 

where the Wilson term and mass term have opposite relative sign compared to the conven- 
tional one. D^g, is the five-dimensional analog of the four dimensional Wilson hopping term 
with 7^ replaced by 75 and U^^x, s) = 1. 

Ds,s' = ^ {(1 - l5)Ss+l,s' + (1 + 75)^5-1,5' - 26s,s'} 

^{(1 -75)5s,L,-l5o,^' + (1 +75)5s,0^L,-l,s'} . (3) 

Note that the five dimensional fermions \E'(x, s) are coupled only to four- dimensional gauge 
fields. The boundaries s = and s = Lg — 1 are anti-periodic and coupled with a weight 

For an appropriate choice of M5, this action has two chiral zero- modes of opposite hand- 
edness, one localized on each boundary of the fifth dimension. To simulate low energy QCD, 
four dimensional quarks are interpolated from the chiral modes, 

q{x) = Pl^{x, 0) + Pr^{x, - 1) , (4) 
q{x) = ^(x, - 1)Pl + ^{x, 0)Pr , (5) 

where Pr l = (l±75)/2. With this definition, one can see that the mj terms in _D^^, directly 
yield the usual four dimensional explicit chiral symmetry breaking term, ■mfq{x)q{x), on the 
four dimensional layers s = and s = Lg — 1. The above definition of the quark fields is the 
simplest one but is not unique. 

In the free theory, for the choice < M5 < 2, the domain wall fermion action corresponds 
to one flavor, and the four- dimensional light quark mass is given as = mfM^{2 — M5) in 
the limit Lg ^ 00 and rrif 0. This situation is approximately unchanged in lattice QCD 
simulations if M5 is simply shifted, M5 — ^ M5 — Mc [T^, where a simple estimate of Mc is 



given in terms of the critical hopping parameter for four-dimensional Wilson fermions 



In the case of domain wall fermions, the SU{Nf) axial transformation can be defined 
vectorially as 

[Q%^{x,s)] =+te{s)X''^{x,s) , (6) 
[Q%^ix,s)] =-teis)^{x,s)X^, (7) 

where s) is a five-dimensional zero mode. Opposite axial charges are assigned to 
fermions in the two half-spaces, so e(s) = 1 for < s < Ls/2 and e(s) = —1 for 
Ls/2 < s < Lg — 1. The reason for this is clear: left and right handed modes are glob- 
ally separated in the extra dimension. 

With this axial transformation on the five dimensional fields, the four dimensional quark 
fields obey the familiar axial transformation: 

[Q%q{x)] =+z^,X-q{x) , (8) 
[Q%q{x)] =-^g(x)75A". (9) 

In general the domain wall fermion action is not invariant under the transformations 
and (|^ in the limit rrif = because of a non- vanishing divergence of the axial current on 
the intermediate layers s = Ls/2 — 1 and s = Ls/2 which gives rise to an extra pseudo- 
scalar density in the axial Ward-Takahashi identity ^j. However, such an anomalous term, 
which provides a chiral symmetry breaking effect due to mixing of the left- and right-handed 
modes, vanishes as — > oo p|. 

Theoretically, this residual breaking effect can be described by an additive quark mass 
m-res in the four- dimensional low energy effective Lagrangian for QCD [|12|, and recent simula- 



tions, in which mres was determined in several ways, appear to confirm this [0. Simulations 
also show that for sufficiently small lattice spacing, mres — as L,, — oo, and at the very 
least, mres is small and accurately known for a wide range of Ls [|12|,|13|. Thus, in this study 



we will henceforth ignore these small effects and assume the chiral limit is = instead 
of rrif = —mres- Finally, we note that the non-perturbative origin of these effects and their 
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relation to certain non-perturbative gauge field configurations is a very interesting and active 
area of research |32 . 



III. BARYON SPECTRUM 
A. Baryon two-point correlator 

The mass of the low-lying baryon B can be extracted from the two-point correla- 
tion function composed of the baryon interpolating operator Ob, which has the appropriate 
quantum numbers specified by the desired state. Let us consider the vacuum expectation 
value of the time-ordered product of interpolating operators. The Euclidean time correlation 
function is projected out at zero spatial momentum through the sum over x: 

GoAt) = E(0|^{^i?(^'^)^B(0,0)}|0) , (10) 

X 

which is dominated by the contribution of the lowest mass state for large Euclidean time, 
Gog{t) ~ exp(— M^ltl). For a spin one-half baryon, of course, the correlation function has 
non-trivial Dirac structure which may be expressed in the form 

Go At) (1 + 74)^(t)^se-^-* + (1 - 74)^(-t)A5e+*^-* , (11) 

large t 

where the positive definite As is proportional to the square of the coupling strength between 
the interpolating operator Ob and the lowest mass state. The first and second terms are 
particle and anti-particle contributions, respectively. It is easy to get only the particle 
contribution by taking the trace with the projection operator P+ = (1 + 74)/2, which we 
abbreviate as ((C>B(t)(5ij(0))) = Tr{P+Go^(t)} hereafter. 

B. Interpolating operators 

Let us focus on the nucleon channel specified by the spin one-half iso-doublet non- 
strange baryons. In order to project out the desired channels, we have to construct interpo- 
lating operators from the quark fields with the appropriate quantum numbers {J = 1/2 and 
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I = 1/2). However, there is considerable freedom in choosing the specific form of the com- 
posite operators. Indeed, there are two possible interpolating operators for the = 1/2"*" 
state, even if we restrict them to contain no derivatives and to belong to the (|, 0) © (0, |) 



chiral multiplet under SU{2)l ® SU{2)r pip 



B^{x) = eabc[ul{x)C-f5db{x)]uc{x) , (12) 
B^{x) = eabc[Uaix)Cdb{x)]-f5Ucix) , (13) 

where abc and ud have usual meanings as color and flavor indices. C is the charge conjugation 
matrix and the superscript T denotes transpose. Here, Dirac indices have been suppressed. 
In Eq.(^2]) and (|13]), the superscript "+" refers to positive parity since these operators trans- 
form as VBi 2{x, t)V^ = +74i?^2(~^) t) under parity. To be precise, the linear combinations 
Bt ± B^ belong to distinct (|, 0) © (0, |) chiral multiplets under SU{2)l © SU{2)r il,|4 



The operator Bi alone is usually used in lattice QCD calculations to extract the nucleon 
ground state since B2 couples only weakly to the ground state due to its vanishing in the 
non-relativistic limit |^. In fact, nobody has succeeded in extracting the nucleon mass 



spectrum from the B2 operator. In our calculation, we also confirm that the ground state 
cannot be extracted from ((-B2^(t)i?2^(0))); however, we have had some success with respect 
to the excited-state mass spectrum which we discuss in Sec IV-C. 

Multiplying the left hand side of the previous positive-parity operators by 75, we obtain 



the interpolating operators with negative parity, J = 1/2 ||19|| : 

B^{x) = 'j^Biix) = €abc[ulix)C'j5dbix)]'j5Ucix) , (14) 
B^ix) = -i^B^{x) = eabc[Ua{x)Cdb{x)]uc{x) , (15) 

since VBi2{x^t)V'^ = — 74i?f2(^3f, t). The important point to notice is the relation between 
the correlation functions of opposite parities, 

GB+it) = --f,GB-{tH, (16) 

since B12 = lf>Bi^2- Eq.(|T^) means that the two-point correlation function of the spin one- 
half baryon can couple to both positive- and negative-parity states. Thus, the general form 
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of the two point function is p5 



Gs-^it) = (1 + j,)9{t)AB^e-^'^^' + (1 - 74)e(-t)AB+e+^B+* 

- (1 + iM-t)AB-e+^'--' - (1 - ^,)e{t)AB-e-''B-i . (17) 

Note that the backward propagating contributions correspond to the anti-particles of the 
forward propagating states with opposite parity. The desired state is obtained by choosing 
the appropriate projection operator, 1 ± 74, and direction of propagation. 

Next, consider a lattice with finite extent T in the time direction and (anti-)periodic 
boundary conditions. Eq. ^ is replaced by 

Gb+W = (1 + iMb+c-''-^' ± (1 - 74)AB+e^*^B+(^-*) 

T (1 + l,)AB^e-''B-(T't) _ (1 _ ^^^^^_^-M,.t ^ (18) 

where the (lower) upper sign stands for (anti-)periodic boundary condition on < t < T. 
The anti-particle of the opposite-parity state can propagate through the time direction 
boundaries, so one faces unwanted contamination from the opposite-parity state in extracting 
the mass of the desired state. This contamination is not unavoidable if a double exponential 
fit is used. However, such fits require very high statistics. This is not a serious issue in the 
measurement of the nucleon ground state since the contamination from the negative-parity 
nucleon is expected to be negligible due to the large mass splitting Mb- — Mb+ ■ Of course, 
this same splitting does affect the extraction of Mb- ■ The problem is resolved by choosing 
appropriate boundary conditions in the time direction to prevent the wrap-around effect, or 
by increasing the time extent T sufficiently and placing the interpolating operators far from 
the boundary. In this study we take the former approach. It is common to employ Dirichlet 
boundary conditions where link valuables in the time direction at t = and t = T — 1 are set 
to zero when calculating quark propagators (for example, see [0]). However, we use a linear 
combination of two quark propagators with periodic and anti-periodic boundary conditions 
in the time direction to produce forward propagating states. Although our approach requires 
two times as many fermion matrix inversions to calculate one quark propagator, it does 
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not suffer from (unknown) reffection effects at the time boundaries induced by Dirichlet 
boundary conditions. 



C. Chiral symmetry and parity doubling 

At the end of this section, let us briefly review how unbroken chiral symmetry imposes 
parity doubling in the hadron spectrum [|14|,^. We will generalize the following argument 
in Appendix A. For the sake of simplicity, we consider a particular transformation of the 
SU{2) chiral symmetry, [Qa,m] = +^75"^ and [Qyi,c?] = —ij^d. Then, one can easily find 
that the two-point correlators Bf and transform as 

'QA,Bt^,ix)Bt^,iO)] =t{l5,Bt^,{x)Bt^,{0)} (19) 

in the chiral limit. Now suppose that the vacuum possesses chiral symmetry: Qa|0) = 0. 
According to Eq.(19), the two-point correlation functions anti-commute with 75, 



{75,G'B±(t)} = 0. (20) 
Immediately, with the help of Eq. ([161) , one finds 

GB+{t)=GBAt) , (21) 



which means that parity doubling arises in the nucleon channel due to chiral symmetry p3 



Of course, in the real world chiral symmetry is spontaneously broken, Qyi|0) 7^ 0, so that such 



parity doubling does not occur in the actual spectrum |T^. In this sense, the spontaneous 
breaking of chiral symmetry is responsible for the absence of parity doubling. Thus, it 
seems important to properly handle chiral symmetry and its spontaneous breaking in order 
to calculate precisely the mass splitting between the nucleon and its parity partner. Finally, 
it is important to note that above argument ignores possible consequences to the 't Hooft 
anomaly condition [p6ip7 |. 
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IV. NUMERICAL RESULTS 
A. Computational details 



We generate quenched QCD configurations on a 16'^ x 32 lattice with the standard 
single-plaquette Wilson action ai [3 = Q/g^ = 6.0. The quark propagator is calculated 
using domain wall fermions with a fifth dimension of Ls=16 sites and a domain wall height 
M5=1.8. Additional details of the simulation can be found in |[12|J38| . Quenched (3 = 6.0 
corresponds to a lattice cut-off of ~ 1.9 GeV from aMp=0. 400(8) in the chiral limit and 



spatial size La ~ 1.7 fm |[1^J38 



We work in Coulomb gauge and calculate quark propagators using wall sources and local 
sinks. We expect that the use of wall sources provides better overlap with the desired states. 
To extract the state with desired parity in the spin one-half baryon spectrum, we construct 
forward (backward) propagating quarks by taking the appropriate linear combination of 
propagators with periodic and anti-periodic boundary conditions in the time direction, as 
mentioned before. This procedure eliminates the backward (forward) propagating opposite 
parity state which wraps around the time boundary. 

In addition, we use two sources for quark propagators on each configuration (tsrc = 5 
and tgj.^ = 27) to increase statistics. After the appropriate parity projections, the correlation 
functions of the corresponding nucleon state can be folded together (averaged as a function 
of distance from the respective source). Here, time-slices are labeled from to 31. To be 
clear, we give the form of the nucleon two-point function with arbitrary source location tgrc 
after eliminating the unwanted contributions across time boundaries: 



Gb+ (t — ts 



;i + 74)AB+e-^^s+(*-w) _ (1 _ ^^)^^_g-M^_(t-w) (T > t > t,,,) , 

;i - 74)^B+e~*'^fl+(*=^^"*) - (1 + 74) ^B-e-^^fl- (tsrc > t > 0) , 



(22) 



which is constructed with the positive-parity interpolating operator (either Bf or -82") • In 
the time range T > t > t^rc, positive- and negative-parity nucleon states are propagating only 
forward in time. This means that there are only particle contributions in this region. Thus, 
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P± = (1 ±74)72 project out positive- and negative-parity states, respectively. Alternatively, 
for tsrc > t > only the anti-particles contribute to the correlation function. In this case, 
P± act in the opposite way to the former case. In either case, all we have to do is calculate 
either Gb+ or Gb- to extract the masses for both positive- and negative-parity states since 
we already know TT{P±GB+(t)} = —Ti{PzpGB-{t)} with the help of Eq.lTBl In fact, we 



verified this relation on each configuration and then used G^+ and G^- to extract masses. 

In our analysis, we use correlation functions in the ranges t > 5 and t < 27, for the 
sources tgrc = 5 and t'^j.^ = 27, respectively. After folding the propagators together as 
described above, we use a single offset from the source, t, defined in the range < t < 27 

We use 405 independent gauge configurations for the lightest two quark masses, mj = 
0.02 and 0.03, 305 configurations for the intermediate ones, rrif = 0.04 and 0.05, and 105 
configurations for the heavier ones, fnf = 0.075 — 0.125. These bare quark masses correspond 
to mass ratios Mt^/Mp 0.59 — 0.90 as shown in Table |. In this calculation, S'f/(2)-isospin 
symmetry is enforced by equating rrif = 171^' = m^p , so that the flavor index will not be 
explicitly displayed hereafter. All calculations were done on the 600 Gfiops QCDSP machine 
at the RIKEN BNL Research Center. 

B. Parity partner of nucleon: A^* 

We first calculate the effective masses Mefj to find appropriate time ranges for fitting. 
The effective mass is defined by 

Meff(t) = ln{C(t)/C(t + l)}, (23) 

where C{t) stands for {{B^2i't) 3^2(0))) ■ We look for a plateau, or time independent region, 
in this quantity to extract the ground state mass. For example. Figures l(a)-l(c) show 
effective masses for the nucleon (-Bi") and its parity partner (-81^2) at mj = 0.03, 0.05, and 
0.10. In Figure 1, the effective mass plot shows a clear plateau for the A^, and one that is not 
as good for the heavier A^*. Statistical uncertainties in Figure 1 are estimated by a single 
elimination jack-knife method. The effective masses for the A^* from both ((i?f(t)i?f (0))) 
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and ((-B^(t)S^(O))) agree well with each other for all three quark masses, except for ruf = 
0.03 and t ^ 6 where there is a small difference outside of statistical errors. 

In Figure 1 (a), the effective mass for the A^* becomes so noisy after t = 9 that the 
value of the two-point correlator is consistent with zero within large errors, so we have left 
these points off of the plot. In addition, the effective mass shows a steep rise with increasing 
uncertainty as t increases. This rise in the A^* effective mass weakens at relatively heavy 
quark masses where the signal becomes stable over 12 time-slices. 

Next we present mass estimates of the and A^* obtained from covariant single exponen- 
tial fits to the corresponding correlators. We fit each correlator from some minimum time- 
slice, tmin, to an appropriate maximum time-slice (tmax = 20 for the A^ and tmax = 10 — 15 
for the A^*). tmax is roughly fixed with reference to the effective mass calculation. To keep 
fitting ranges as wide as possible, tmin is reduced from tmax — 2 until x^/^df > 1-5 where 
A^DF denotes the degrees of freedom in the fit. Fitting details are given in Tables |Il-|V[ All 
of our fits have confidence-level larger than 0.2 and estimates from the weighted average of 
the effective mass agree with the fitted masses within errors. A summary of our A^ and A^* 
masses is given in Tables [ni-[rV|. 

We can roughly estimate the systematic error coming from the choice of fitting range 
by varying tmin- Let us determine the difference between our final fits and fits where the 
smallest time slice is not included. In the case of (the nucleon), resulting errors are much 
smaller than the statistical errors of the final fits. On the other hand, for both i?f and B2 , 
the estimated systematic errors are comparable to the statistical errors, while for at 
TTLf = 0.04 and 0.05 the systematic error is a factor of two larger than the statistical one. 
This uncertainty is directly related to the rise of the A^* effective mass, or poor plateau, 
in those cases where the correlation function was not statistically well resolved at larger 
times. This, in turn, has forced us to extract masses from the smaller t region which may 
be contaminated by excited states. 

In Figure 2 we show the low-lying nucleon spectrum as a function of the quark mass, rrif. 
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The nucleoli mass is extracted from Bf. We omit the point at mj = 0.02 for the operator 
B2 since a good plateau in the effective mass plot is absent. The A^* mass estimates from 
((-Bf (t)-Bf (0))) and {{B2 {t) B2 (0))) agree with each other within errors in the whole quark 



mass range, as expected from their common quantum numbers |T^. Note that this result 
disagrees with that obtained in : we find no discrepancy between masses extracted from 
((fif (t)i?r(0))) and {{B2{t)B2{0))). However, recent resuhs in ^ and ^\ are m good 
agreement with ours. We also obtain the same mass for the A^* from a mixed correlator 

{{B^it)B2iO) + B2{t)B^m)- 



The most remarkable feature in Figure 2, which was first reported in [|15|, is that the 
N-N* mass splitting is observed over the whole range of quark mass values and grows as 
the quark mass is decreased. To illustrate this point clearly, we compare two mass ratios in 
Figure 3, one from the baryon parity partners Mi\f*/Mj\f and the other from pseudo-scalar 
and vector mesons M-^/Mp. Experimental points are marked with stars, corresponding 
to non-strange (left) and strange (right) sectors. In the strange sector we use S"*" and S(1750) 
as baryon parity partners and K and K* for the mesons. The baryon mass ratio clearly 
grows with decreasing meson mass ratio, toward the experimental values [0. We did not 
include the charm sector [M^i/Mj:,* ~ 0.93) since the parity partner of S++(2455) is not 
measured experimentally. On the other hand, from our results we estimate the mass of this 
state to be roughly 2.7 GeV (we have used a simple linear in the quark mass ansatz to 
extract the mass from our degenerate quark data). 

Finally, we evaluate the N and A^* masses in the chiral limit. Taking a simple linear 
extrapolation in the four lightest quark masses for B^ and B^ , we find MAr=0.57(l) and 
M7v*=0.85(5) in lattice units. Setting the scale from the calculated p-meson mass |jl2|, we 



obtain ^1.1 GeV and M^* ^ 1.6 GeV in the chiral limit. If we use the scale set 
by the calculated nucleon mass, we obtain Mjy* ~ 1.5 GeV. Either way the A^* mass is 
in good agreement with the experimental values within about 5-10 %. The above errors 
do not include systematic uncertainties due to finite volume, non-zero lattice spacing, and 



16 



quenching effects. Studies of such systematic errors will be addressed in future calculations. 

Recently it was suggested that the A^* propagator in quenched QCD may exhibit 
non-analytic chiral behavior associated with the anomalous contribution of the A^" inter- 
mediate state.0 This suggestion is inspired by recent articles ||41|,^ . However, non-analytic 



effects which arise in the vicinity of the massless quark limit are hard to detect in our data 
since our calculation is not very close to this limit. In addition, such subtle effects could 
easily be mimicked by finite volume effects. The volume used in our study is still rather 
small. Furthermore, the statistical errors in our study are large enough, especially for the 
N*, that complicated fits beyond a naive linear one would not yield meaningful parameters 
or x^- 

C. Unconventional nucleon operator 

As mentioned earlier, the unconventional operator B2 vanishes in the non-relativistic 
limit, so one may expect that it couples with negligible weight to the nucleon ground state 
near this limit. Indeed, no one has succeeded in extracting the ground state mass signal 
from {{B2 (t) B2 (0))) in lattice QCD [ p8| , p^ . Then, we expect that an approximately zero 
overlap on the ground state provides the possibility of direct access to the excited state using 
the B2 operator. This prospect should hold as long as the lattice operator has the same 
symmetries as the continuum one. 

In the massless quark limit, the combinations Bi ± B2 do not mix due to different 
chiral structure. In perturbation theory, each is multiplicatively renormalized with the same 
renormalization factor so that Bi and B2 also do not mix. However, conventional lattice 
fermions give rise to mixing through explicit chiral symmetry breaking Thus, B2 will 



couple to the ground state through unwanted mixing with Bi . On the other hand, the 



^ Recall that in the quenched approximation the ry' is also a pseudo-Goldstone boson. We stress 
that the quenched "r/'-A^" intermediate state would give a negative metric contribution correspond- 



ing to a "quenched chiral loop" artifact [41|. Needless to say, this spurious intermediate state, which 



is a source of unitarity violation, is not associated with the physical decay process [41|. 
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explicit breaking of chiral symmetry in domain wall fermions is highly suppressed, so there 
is hope that we may extract the positive-parity excited state of the nucleon from B2 as the 
asymptotic state. 

Let us first compare the effective mass plots of {{Bi(t)Bi{0))) and ((i?^(t)5^(0))) corre- 
lators (see Figures 4a-4c). The correlator {{B2{t)B2{0))) is considerably noisier so that only 
time-slices near the source are useful. Nevertheless, the effective mass from ((52^(t)i?2^(0))) 
yields a plateau albeit with large statistical errors. The plateau becomes more satisfactory 
for heavier quark mass. These plateaus are obviously different from those extracted from the 
B^ correlator. When we apply the single exponential fit to the two-point correlation func- 
tion composed of B2, we obtain a mass that is quite large compared to the mass extracted 
from {{Bi (t) Bi [0))) . In pO| similar results from the B2 correlator have been found. 



As mentioned earlier, an explanation of the above result [jT5| is that B2 has negligible 
overlap with the nucleon ground state since it does not have a non-relativistic limit and thus 
provides a direct signal for the positive-parity excited state of the nucleon. Of course, this 
explanation is valid only in the heavy valence-quark mass limit. 

We investigate further the possibility that (0|-B^|A^) ~ through the calculation of the 
mixed correlation function {{Bi(t)B2{0) + i?2^(t)i?j^(0))). If it is true that the overlap with 
the nucleon (0|-B^|A^) becomes small with increasing valence-quark mass, mass estimates 
from the mixed correlator {{Bi(t)B2{0) + B2{t)Bi{0))) should be consistent with the 
nucleon for lighter quark mass and the positive-parity excited state for heavier quark mass. 
Such behavior is evident in Figure 5. We stress that the single particle fit to the mixed 
correlation function does not in general yield the mass of an asymptotic state, only in the 
limits discussed above, which is clear from the figure. We also note that for improved Wilson 
fermions the mixed correlation function yields the ground state mass for all quark masses. 



even heavy ones Thus it appears in that case that there is still significant mixing of 
B^ and B^. 

As a result, it is possible to identify B2 with the positive-parity excited nucleon (A^') for 
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heavy quarks {e.g. nif > 0.075 in our study) pf. Indeed, we see |(0|S^|A^)/(0|S^|A^')P < 
10^'^ from double exponential fits for ((i?2'(t)-B2^(0))) at nif = 0.10 and 0.125. Of course, 
this feature weakens in the lighter quark mass region (from around rrif = 0.05). However, 
even for rrif < 0.075, mass estimates from B2 are still considerably larger than the nucleon 
mass, so we cannot rule out the possibility that B2 provides a signal for the positive-parity 
excited nucleon at even lighter quark mass values. 

D. Diagonalization method for 2x2 matrix correlator 

Let us consider the N' in more detail. The mass of the excited nucleon may be 
obtained from a two state fit to the Bf correlator which, however, requires large statistics 
and that neither the ground state or the first excited state lies close to any other state. 
Attempts to extract the A^' mass using this method have failed to reproduce the observed 
mass ||2^. In our fits the Hessian matrix often becomes singular, so we cannot extract the 
excited state mass from double exponential fits. Instead, we take an alternative approach 
proposed in |^^. First, we define the 2x2 matrix correlator C{t) using the two distinct 
baryon operators 



C{t) 



Cll(t) Ci2(t) 
C2l(t) C22(t) 



(24) 



where Cij{t) = {{Bf{t)B^{0))). Next we write C{t) in terms of a transfer matrix A(t,to), 

C{t)ij = X{t,to)C{to)^ , (25) 

where to is fixed and t > to. If only two states are propagating in a given system, the masses 
of the two states {Ei > Eq) are given by the eigenvalues of the transfer matrix: 

A„(t,to)=e-(*-*«)^'^ (« = 0,1), (26) 

where E^ is independent of to- The smaller eigenvalue (Ai) and larger eigenvalue (Ao) refer to 
the masses of the excited state and the ground state respectively. In general, the system may 
have more than two states. Thus, we assume that two states become effectively dominant 
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for an appropriately large time-slice to, which can be determined by checking the sensitivity 
of Ea with respect to variations of to- 

We calculate the eigenvalues Xait,to) of C(to)~^C(t) for to = 3, and then evaluate the 
effective masses of and A^' from the larger and smaller eigenvalues, respectively, with 
the statistical errors coming from the jack-knife method. Using these error estimates, we 
determine the masses in Table VI from weighted averages over the time-slice ranges listed 
there. Figure 6 shows that the results are quite consistent with the masses determined from 
single exponential fits to ((S^ (t)S^ (0))) and {{Bi{t)B^{0))). We have checked that the 
variations of to around to = 3 do not significantly affect the effective masses. For rrif < 0.03, 
we unfortunately could not extract the mass of the first excited nucleon since a good plateau 
in the effective mass plot is absent. We note that for mj = 0.04 there is still a large splitting 
in the eigenvalues, indicating the overlap of B2 with the nucleon remains small. 

In contrast to the positive-parity state, the eigenvalues in the negative-parity case appear 
degenerate, as shown in Figure 7. This will be discussed in more detail in the next section. 

As we have seen, the first excited state may die out so quickly that only a few time-slices 
are available for evaluating the mass, even if the excited state is well separated from the 
ground state. To circumvent this quick damping, simulations performed on an anisotropic 
lattice where the temporal lattice spacing is finer than the spatial one may be useful. In 
fact, a recent lattice study [^] has shown this to be an effective way to extract masses of 
nucleon excited states. 

E. Comparison with experiment 

In the physical spectrum, we have another negative-parity state, A^*(1650), which is 
just above the lowest state A^*(1535). Although these states are quite close to each other, 
they are easily distinguished due to a peculiar decay mode of A^*(1535). It is well known that 
the decay rate of A^*(1535) —>■ N + ri is comparable with A^*(1535) —>■ N + 7t even though the 
Nn decay mode is kinematically favored over Nrj. For the case of A^*(1650) the A^vr decay 
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is dominant; the branching ratio of Nrj is only a few percent. Thus, A^*(1535) N + % 
seems to be anomalously suppressed. However, the corresponding states cannot be easily 
distinguished without knowledge of decay patterns if they lie close to each other. 

As mentioned, for the negative-parity case our data shows that both eigenvalues of the 
transfer matrix are the same within statistical errors (see Figure 7). Of course, this does not 
rule out the possibility that the splitting might become clear with more statistics and in the 
lighter quark mass region. That is because the central value of the difference between the 
two eigenvalues increases in the lighter quark mass region. Furthermore, the value of about 
100 MeV at irif = 0.03 seems to be large enough to reproduce experimental splitting in the 
chiral limit. 

Nevertheless, unlike the positive-parity case, we have no indication of the presence of 
two independent negative-parity states in our calculation where the mixed correlator yields 
a consistent mass with both and B2 interpolating operators, as listed in Table VI. Thus, 
our combined results from all data for the negative-parity nucleon allows two other possibil- 
ities regarding A^*(1650). One is that A^*(1535) and A^*(1650) are completely degenerate or 
not independent in the quenched calculation. This situation resembles naive quark models 
and the MIT bag model where A^*(1535) and A^*(1650) are degenerate if we neglect the spin 
dependent interaction. In this case, the analysis through the diagonalization of the 2x2 
correlator is no longer helpful. 

A Second possibility is that A^*(1650) is simply missing in our calculation, i.e. both 
Bi and B2 operators may not couple, or couple weakly, to the A^*(1650) state. This may 
be related to an argument in |Q where the authors show that the desired vanishing of 
the phenomenological ttNN* coupling for the A^*(1535) state results simply from the chiral 
transformation properties of particular interpolating operators i?f and B2 with the help of 
the soft pion limit. If their argument is relevant to the suppression of the A^vr decay, it is 
possible that -Bf and B2 couple strongly to A^*(1535) but not to A^*(1650). However, their 
argument is no longer applicable to a third interpolating operator 
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-^3 M = {u^C(^al3u)(TapJ^d - {u^ C (Japd)(Jap'^ , (27) 



which belongs to the chiral multiplet (|, 1) © (1, |) and has no derivative ||33| . |3^ . The By, ^ 



couples to both J = 3/2 and J = 1/2 states. In this case, A^*(1650) might be extracted 
from the -83"^ operator since there is no reason for lack of coupling to the A^*(1650) state. 

Finally, we would like to mention a remaining puzzle. We have done the first successful 
lattice calculation of both N* and A^' spectra. As for the N' , we have reliable data only 
for relatively large values of mj. Comparing the A^' mass with the A^* mass, we find that 
the ordering of A^' and A^* is inverted compared to experiment. Furthermore, the level 
spacing between N-N* and N*-N' is almost even ||T5[. What we see here closely resembles 



the wrong ordering problem of the excited nucleon spectrum in naive quark models and the 
MIT bag model. However, if this result is true for heavy quarks, we can make an important 
prediction: the first excited state of the spin one-half Sc is a negative-parity state rather 
than a positive-parity state. 

Unfortunately, our statistics did not allow the computation of the A^' mass in the light- 
quark region {rrif < 0.03). In addition, our data show no evidence for the possibility of 
switching the level ordering between A^* and A^' towards the chiral limit. Indeed, we did not 
observe that the A^' mass decreased faster than the A^* mass with decreasing quark mass. 
However, finite volume effects for higher excited states may be more serious than for lower 
excited states. According to naive quark models or the MIT bag model, the A^' state is 
radially excited in contrast to the nucleon ground state and the A^* state. We note that the 
trend in the three heaviest quark mass points is at least consistent with the switching of the 
ordering in the chiral limit. Needless to say, the quenched approximation may also play a 
role in this puzzle. In summary, our results do not rule out the possibility of switching the 
ordering between A^* and A^' near the chiral limit, and to solve this remaining puzzle we 
need further systematic calculations toward the chiral limit. 
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V. CONCLUSION 



We have studied the mass spectrum of the parity partner of the nucleon in quenched 
lattice QCD using domain wall fermions which preserve chiral symmetry to a high degree 
in lattice simulations. Most importantly we demonstrated that this method is capable of 
calculating the mass of the negative-parity A^* state in the spin one-half and isospin one-half 
baryon sector. 

We made a systematic investigation of the A^* spectrum by using two distinct interpo- 
lating operators, and B2. We found the N* mass estimates extracted from them agree 
with each other. In practice the -Bf correlator is more convenient than B2 in extracting the 
A^* mass because it is less noisy, especially in the light quark mass region. 

We have found a definite mass splitting between N and A^* states in the whole quark mass 
range we studied. Furthermore, this splitting grows with decreasing quark mass. The A^* 
mass and the N-N* mass splitting in the chiral limit obtained by extrapolation are consistent 
with the experimental value within about 5-10%, depending on the mass we choose to set 
the lattice spacing, the nucleon or the p meson. These results have been confirmed by 



subsequent lattice calculations using improved Wilson fermions p0|j21[] . Needless to say this 
is very encouraging for further investigations of A^* physics using lattice QCD simulations. 

In contrast to the negative parity operators, the positive parity operators, B^ and B2, 
yield distinct mass signals. From the B^ operator we obtained a clean signal for the nucleon 
ground state, while from B2 we always found a heavier mass. This is probably because 
the latter vanishes in the non-relativistic limit and thus has small overlap with the ground 
state. Indeed, we confirmed numerically that (0|i?2^|A^) ^ and that B2 yields the mass 
signal of the positive-parity excited state in the heavy quark region. We found that in 
the heavy quark mass region the mass of this excited state is about twice as high above 
the ground state mass as the A^* mass. This property of the B2 operator was confirmed by 
comparing the results with the excited-state spectrum obtained from the diagonalization of a 
2x2 correlation matrix constructed from both interpolating operators and also by examining 
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single particle fits to the mixed correlation function. 

We did not resolve a long-standing puzzle regarding the excited-state spectrum of the 
nucleon, the inversion of the positive- and negative- parity first excited states with respect 
to experiment. However, we beheve that the level switching between N* and N' must occur 
close to the chiral limit, though there is scant evidence in our results that this might happen. 
Taking our calculation at face value leads to the prediction that the first excited state of Ec 
has = 1/2-. 

Needless to say, more work needs to be accomplished in order to achieve a fully systematic 
calculation of the excited nucleon spectrum. In the near future we plan to increase our 
volume and statistics to further explore the chiral limit, and to check that lattice spacing 
errors are small by running at several gauge couphngs. We expect that the latter effect is 
small due to the chiral symmetry properties of domain wall fermions. A long term goal is 
to include dynamical fermions. 
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APPENDIX A: GENERAL CHIRAL TRANSFORMATION AND PARITY 

PARTNERS 



To extend our discussion of Sec.III-C to general chiral transformations, we prepare 
iso-doublet operators for Bi and B2, 

' ^ ' \ eabc[dl{x)C^,Ub{x)]d,{x) ) ' ^^""^ 

„+f . _ ( eabc[ul{x)Cdb{x)]-i^Uc{x) \ , . 

\ eabc\dl{x)Cub{x)\^^dlx) ) ■ ^^^^ 

The upper component corresponds to the proton and the lower component corresponds to 
the neutron. Also we can define iso-doublet operators for the negative parity nucleons as 
i3]^2 = 75-^1^2) respectively. Here we consider a general transformation of SU (2)v and SU (2)^ 
symmetry, 

\Q%,q{x)\ =iT,q{x) , (30) 
\Q\,q{x)\ =iT,^^q{x) , (31) 

where q = {u, d)'^ and is a Pauli matrix. Under these transformations, the two-point 
correlator composed of either Bf or Bf should transform as 

Q^, Bt^,ix)Bt^,iO)] = ^ [ra, Bt^,ix)Bt^M] , (32) 
Q%B^^,{x)Bt,{0)] =^{l5r.,BUx)BU0)} . (33) 

Eq.(^3]) along with the fact that Qy|0) = tells us that the nucleon two-point correlation 
function should be commute with r^, 

[Ta,Gis±{t)]=0, (34) 

which means that G'g+ (t) has only diagonal elements of iso-spin indices. In other words, the 
proton state and the neutron state are eigenstates of iso-spin, of course. If Q^|0) = 0, it 
turns out that 



{l5ra,Gs±{t)} = r„{75,Ge±(t)} = 
25 



(35) 



Thus, we can obtain the strict relation in terms of two-point correlators between the nucleon 
and its parity partner 

Gs+it) ^ Gs-{t) (36) 
in the non-broken phase of chiral symmetry. 
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FIGURE CAPTIONS 



FIG.IA The effective mass of the nucleon from the {{B^{t)B^{0))) correlator (x) 

and its parity partner from ((Bf (0))) (□) and {{B2 it)B2 (0))) {O) 
correlators on an ensemble of 405 configurations at rrif = 0.03. The solid 
lines and dashed lines represent each fitted mass and its statistical error. 



FIG. IB The effective mass of the nucleon from the {{Bi{t)Bi{0))) correlator (x) 

and its parity partner from ((5f (t)5f (0))) (□) and {{B2 (t) B^ (0))) {O) 
correlators on an ensemble of 305 configurations at = 0.05. The solid 
lines and dashed lines represent each fitted mass and its statistical error. 



FIG.IC The effective mass of the nucleon from the {{Bi(t)Bi{0))) correlator (x) 
and its parity partner from (t)5f (0))) (□) and {{B2 {t)B2 (0))) {<>) 
correlators on an ensemble of 105 configurations at rrif = 0.10. The solid 
lines and dashed lines represent each fitted mass and its statistical error. 

FIG. 2. iV (x) and A^* (□ and O) masses versus the quark mass in lattice 
units (a^^ ^ 1.9GeV from aMp=0. 400(8) in the chiral limit). The corre- 
sponding experimental values for N and N* are marked with lower and 
upper stars. The N-N* mass splitting is clearly observed. Symbols (x, 
□ and O) are defined as in Figure 1. 

FIG. 3. Mass ratio of the negative-parity excited-state and positive-parity 
ground-state baryons versus mass ratio of the pseudoscalar meson 
and vector meson. All calculations are done for three degenerate 
valence-quarks. Ratios are calculated using fitted masses from the 
((Sf (i)Sr(O))) (□) and ((S2~(i)S2~(0))) (O) correlators. Experimental 
points are denoted by stars. 
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FIG.4A The effective mass from the ((Sj^(i)Si+(0))) (x) and {{B}{t)Bt{0))) (o) 
correlators on an ensemble of 405 configurations at m/ = 0.03. The solid 
lines and dashed lines represent each fitted mass and its statistical error. 

FIG.4B The effective mass from the {{Bt{t)Bt(0))) (x) and {{Bt(t)Bt(0))) (o) 
correlators on an ensemble of 305 configurations ainif — 0.05. The sohd 
lines and dashed lines represent each fitted mass and its statistical error. 

FIG.4C The effective mass from the {{Bt(t)BtiO))) (x) and {{Bt{t)Bt{0))) (o) 
correlators on an ensemble of 105 configurations at m/ = 0.10. The solid 
lines and dashed hues represent each fitted mass and its statistical error. 

FIG. 5. The fitted mass from {{Bt(t)Bt{0))) (x), {{Bt(t)Bt{0))) (o) and mixed 
type {{Bi{t)B^{0) + B:^{t)Bt{0))) {O) correlators. The corresponding 
experimental values for and N' are marked with lower and upper 
stars. Note, the values extracted from the mixed correlation function do 
not represent the mass of an actual asymptotic state, except possibly in 
the light and heavy quark limits. 

FIG. 6. Comparison of the fitted mass from ((-B^(t)i?^(0))) (o) and the esti- 
mated mass from the average effective mass of the smaller eigenvalue 
of the transfer matrix (•). The symbol x corresponds to the nu- 
cleon ground-state mass evaluated from the larger eigenvalue of the 
transfer matrix, which is quite consistent with the fitted mass from 

{{Btmtm). 
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FIG. 7. Symbols □ and O correspond to the the estimated mass from the average 
effective mass of the larger and smaller eigenvalues of the transfer ma- 
trix for the negative parity state. The eigenvalues are degenerate within 
errors. Symbols x and o are defined as in Figure 6. The correspond- 
ing experimental values for 7V(939), A'"'(1440) and A^*(1535) are marked 
with lower, middle and upper stars. The ordering of the negative-parity 
nucleon (iV*) and the positive-parity excited nucleon (A^') is inverted 
relative to experiment. 
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TABLES 











MjMp 


jj^ of configs. 


0.020 


0.2687(24) 


0.4530(62) 


0.593(14) 


98 


0.030 


0.3224(21) 


0.4814(45) 


0.670(11) 


98 


0.040 


0.3691(19) 


0.5126(42) 


0.720(10) 


98 


0.050 


0.4116(18) 


0.5395(36) 


0.763(8) 


98 


0.075 


0.5080(17) 


0.6088(38) 


0.834(8) 


98 


0.100 


0.5962(18) 


0.6797(36) 


0.877(7) 


98 


0.125 


0.6774(17) 


0.7483(33) 


0.905(6) 


98 



TABLE I. Single exponential fit of tt and p two-point correlators at /? = 6.0 on a 16^ X 32 X 16 
lattice with M5 = 1.8. The masses for anif < 0.05 agree within errors with those given in 
Ref. [12 1 which were calculated from different correlation functions on a slightly different set of 
configurations. 



state (J^) 


anif 


aMAr(stat.)(sys.) 


xVAdf 


^min 


^max 


conf. level 


# of configs. 


N (1/2+) 


0.020 


0.654(12)(14) 


1.00 


12 


20 


0.43 


405 




0.030 


0.716(5)(0) 


1.09 


10 


20 


0.37 


405 




0.040 


0.754(6)(1) 


1.08 


11 


20 


0.38 


305 




0.050 


0.805(5)(2) 


1.17 


11 


20 


0.31 


305 




0.075 


0.929(7)(4) 


0.81 


12 


20 


0.58 


105 




0.100 


1.045(5)(3) 


0.94 


12 


20 


0.47 


105 




0.125 


1.162(6)(4) 


0.84 


15 


20 


0.50 


105 



TABLE IL Single exponential fit of the nucleon two-point correlator {{Bf {t)Bf {0))) at (3 = 6.0 
on a 16'^ X 32 X 16 lattice with M5 = 1.8. The systematic error is estimated from the change in 
the fitted mass when tmin is increased by one. The masses for anif < 0.05 agree within errors 
with those given in Ref. |12] which were calculated from different correlation functions on a slightly 
different set of configurations. 
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state (J^) 


aiUf 


aMj\f* (stat.)(sys.) 




^min 


^max 


conf. level 


# of configs. 


N* (1/2-) 


0.020 


0.935(34) (40) 


1.29 


6 


10 


0.28 


405 




0.030 


0.959(24)(24) 


1.16 


7 


12 


0.33 


405 




0.040 


1.018(28)(56) 


1.49 


8 


12 


0.22 


305 




0.050 


1.048(20)(38) 


1.26 


8 


13 


0.28 


305 




0.075 


1.104(13)(18) 


1.27 


7 


13 


0.27 


105 




0.100 


1.203(10)(12) 


0.99 


7 


12 


0.41 


105 




0.125 


1.303(8)(8) 


0.92 


7 


12 


0.45 


105 




TABLE III. Sin^ 


;le exponential fit 


of the nej 


^ative-parity 


nucleon two-point correlator 




(0))) at /? 


= 6.0 on a 16^ x 32 x 16 lattice with M5 = 


1.8. The systematic error is 


estimated from the change in the fitted mass when tmin is increased by one. 






state {J^) 


aruf 


aMjy* (stat.)(sys.) 




^min 


^max 


conf. level 


# of configs. 


N* (1/2-) 


0.030 


0.982(21)(15) 


1.10 


5 


10 


0.35 


405 




0.040 


0.993(16)(17) 


0.81 


5 


10 


0.52 


305 




0.050 


1.029(15)(16) 


1.00 


6 


10 


0.39 


305 




0.075 


1.092(12)(10) 


0.75 


5 


11 


0.59 


105 




0.100 


1.193(11)(16) 


0.82 


6 


12 


0.54 


105 




0.125 


1.308(13)(9) 


0.84 


7 


15 


0.55 


105 




TABLE IV. Sin^ 


^le exponential fit 


of the nef 


jative-parity 


nucleon two-point correlator 


{{B^it)B^{0))) at P 


= 6.0 on a 16^ x 32 x 16 lattice with M5 = 


1.8. The systematic error is 


estimated from the change in the fitted mass when tmin is increased by one. 






state (J^) 


am,f 


aMjv' (stat.)(sys.) 


xVA^df 






conf. level 


# of configs. 


N' (1/2+) 


0.030 


1.170(50)(108) 


1.47 


5 


8 


0.23 


405 




0.040 


1.191(36)(69) 


1.28 


5 


11 


0.27 


305 




0.050 


1.244(43)(148) 


1.49 


6 


11 


0.20 


305 




0.075 


1.303(39) (58) 


0.43 


6 


10 


0.73 


105 




0.100 


1.387(19)(2) 


0.63 


5 


10 


0.64 


105 




0.125 


1.484(16)(6) 


1.00 


5 


11 


0.42 


105 



TABLE V. Single exponential fit of the unconventional nucleon two-point correlator 
{{B^{t)B^{0))) at /? = 6.0 on a 16^ x 32 x 16 lattice with M5 = 1.8. The systematic error is 
estimated from the change in the fitted mass when tmin is increased by one. 
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arrif 


Mixed type 


[^miri) ^max] 


Eo 


[^min) ^max] 


El 


[^rniii) ^max] 


# of configs. 


1/2+ 


0.020 


0.704(67) 


[6,11] 


0.649(10) 


[ 6,13] 


N/A 




200 




0.030 


0.795(57) 


[5,12] 


0.703(9) 


[ 7,15] 


N/A 




200 




0.040 


0.873(56) 


[5,12] 


0.755(9) 


[ 7,20] 


1.264(84) 


[5, 8] 


200 




0.050 


0.942(54) 


[5,11] 


0.808(8) 


[ 8,20] 


1.247(67) 


[5,10] 


200 




0.075 


1.198(115) 


[5,10] 


0.927(11) 


[11,20] 


1.256(39) 


[4,10] 


81 




0.100 


1.405(54) 


[3,10] 


1.048(9) 


[12,20] 


1.382(39) 


[5,10] 


81 




0.125 


1.504(83) 


[4,10] 


1.161(7) 


[13,20] 


1.479(33) 


[5,10] 


81 


1/2- 


0.020 


N/A 




0.901(50) 


[5,10] 


N/A 




200 




0.030 


0.992(47) 


[5,11] 


0.938(29) 


[5,10] 


1.004(49) 


[4,10] 


200 




0.040 


0.979(24) 


[5,11] 


0.996(27) 


[6,10] 


1.021(35) 


[5,10] 


200 




0.050 


1.005(16) 


[5,11] 


1.054(29) 


[7,11] 


1.042(26) 


[5,11] 


200 




0.075 


1.092(14) 


[5,12] 


1.105(18) 


[5,10] 


1.102(22) 


[5,10] 


81 




0.100 


1.185(10) 


[5,10] 


1.197(18) 


[5,12] 


1.203(16) 


[5,12] 


81 




0.125 


1.282(9) 


[5,11] 


1.292(16) 


[5,12] 


1.303(13) 


[5,12] 


81 



TABLE VI. The third column hsts the fitted masses of the mixed type correlators. Fifth and 
seventh columns list the weighted averages of the effective mass derived from the larger eigenvalue 
and the smaller eigenvalue of the transfer matrix Xa{t, to = 3) induced by the 2x2 matrix correlator. 
The fits corresponding to column three all have x^/^df < 1-5 For the positive-parity state, the 
larger and smaller eigenvalue of the transfer matrix are clearly distinguishable. They correspond 
to the nucleon ground state (N) and the first excited nucleon {N') respectively. However, for the 
negative parity state, both eigenvalues are degenerate within errors. 
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FIG. la (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG. lb (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG.lc (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG.2 (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG. 3 (Phys.Rev.D) Shoichi Sasaki et al. 




0.4 0.6 

Mps / Mv 



1.0 



37 



FIG. 4a (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG.4b (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG.4c (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG. 5 (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG. 6 (Phys.Rev.D) Shoichi Sasaki et al. 
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FIG. 7 (Phys.Rev.D) Shoichi Sasaki et al. 
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